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Abstract 

Within Euclidean lattice field theory an exact equivalence between the 
one-flavour 2D Thirring model with Wilson fermions and Wilson parameter 
r = 1 to a two-colour loop model on the square lattice is established. For non- 
interacting fermions this model reduces to an exactly solved loop model which 
is known to be a free fermion model. The two-colour loop model equivalent 
to the Thirring model can also be understood as a 4-state 49-vertex model. 



^E-mail: scharnhOlinde . physik . hu-berlin . de 

* Present address 



1 Introduction 



Four-fermion models on the lattice and in the continuum (in various dimensions) 
have received a great deal of attention in recent years. This interest is related to 
the problem of dynamical symmetry breaking |]J, p] as well as to the possible role 
of four-fermion interactions in QED in 4D ||, [|3J. However, the reliable investiga- 
tion of these models requires the application of non-perturbative and exact methods 
which are hard to come by. Within lattice field theory the investigation of fermionic 
theories is hampered by the notorious difficulties encountered in their numerical 
simulation, in particular, when the number of fermion species involved is odd. The 
development of new methods for this class of theories is therefore highly desirable. 

Following a method developed by Salmhofer || in the exact study of the one- 
flavour strong-coupling lattice Schwinger model with Wilson fermions, we have re- 
cently suggested || that this method can also be applied to a wide class of purely 
fermionic theories with Wilson fermions in various dimensions. The study of the 
one- and two-flavour strong-coupling Schwinger models has lead to exact equiva- 
lences to statistical models which can be understood as either q-state vertex models 
or loop models 0, @. In contrast to strongly coupled gauge theories with stag- 
gered fermions, which lead to pure monomer-dimer systems , the consideration of 
Wilson fermions results in a different class of model, as these examples demonstrate. 



In the present article, as a first step towards the application of the new method 
to four-fermion theories, within Euclidean lattice field theory we investigate the one- 
flavour 2D Thirring model with Wilson fermions. In the past decades the one-flavour 
Thirring model in 2D has been widely studied in the continuum as well as within 
the Hamiltonian lattice approach. The massive Thirring model in the continuum is 
known to be quantum integrable (||, see also ||) and to be equivalent to the sine- 
Gordon model [ITn] . Within the Hamiltonian lattice approach equivalences of the 



Thirring model to other statistical models have been found previously [11], [T2| (see 



also the references cited in ||). These models are related to staggered fermions. The 
only attempt so far to study the one-flavour Thirring model with multi-component 



spinors has been made by Ishida and Saito JT3[. They however consider the case of 
naive fermions ('Wilson fermions with Wilson parameter r = 0') only, which turns 
out to be much more involved than the case of Wilson fermions with r = 1 we are 
going to investigate in the present article. 



The method we will employ is based on the fact that the Grassmann integration 
on the lattice can be performed exactly. The Grassmann integration will consecu- 
tively be carried out on the even and odd sublattices. In this process the possible 
contributions to the partition function will be analyzed (for the method see also || , 
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3] ) . In order to simplify the required algebraic manipulations with Grassmann vari- 



ables we make use of a Mathematica program |14] specially designed for the present 
purpose. Since we consider only one fermion species, throughout the paper we pay 
careful attention to the ordering problem for the Grassmann variables involved in 
order to keep track of any minus signs which possibly might lead to negative con- 
tributions to the partition function. 

The plan of the article is as follows. The main investigation is presented in 
section 2 which is divided into several subsections. This division is introduced in 
order to give the reader easier access to the different steps in establishing the trans- 
formation from the one-flavour lattice Thirring model with Wilson fermions to the 
equivalent statistical model. The structure of section 2 is explained in further detail 
at the end of its introductory part (before subsection 2.1 starts). Section 3 contains 
some discussion and conclusions. Finally, the proof of a mathematical result needed 
in subsection 2.2 is deferred to an Appendix. 



2 The loop model equivalence 

The partition function Z\ of the one-flavour Thirring model with Wilson fermions 
on the square lattice A is given by 



Z A = J DtpDtP e~ s , (1) 

where DipDip = TI^ga c^/M^) dipi(x) dip2{x) d^ix) denotes the multiple Grassmann 
integration on the lattice (The subscripts are spinor indices.). The action S is defined 
by 

- M$(xWx) - % J^J^x)) (2) 



with 



J fi( x ) = ^(xh^(x) , (3) 
J M (x)J M (x) = 4 ^(x^ix^iix^ix) . (4) 



In eq. (^) the Wilson parameter r has been set to 1; the hopping parameter k is 
consequently given by k = 1/2M. Eq. ([!]) can be written as follows. 



J2 I Dtp Dtp exp ( Mip(x)i)(c 

), 1,2,3,4} V xeA 

exp 2G 5^ ^i(x)^ 2 (x)il;i(x)iff2(x) 
\ xeA / 



x II ^y[^(x + e M )T / ( + V(a;) + ^(x)T(-)^(x + e Ai )] fci . (5) 

The occupation numbers fcj cannot exceed 4 because on each lattice site only 4 
Grassmann variables are present. Explicit calculation however shows that by virtue 
of the nilpotency of the Grassmann variables k\ cannot exceed 2 (This is closely 
related to the fact that we have employed r = 1.). As in JJ, we choose 71 = 03, 
72 = ®\ [p% are Pauli matrices.) and the projection operators T^> = (1 + e r y fl )/2, 

e G {-1,1} then read t} +) = diag{l,0}, = diag{0, 1}, T 2 (e) = It is 

now advantageous in order to simplify further results to define combinations of the 
original ip fields in terms of which T^ 1 is diagonal 



1 T 

x ( x ) = Uip(x) = -^(^(x) +ip 2 (x), -01 (x) -V2 (x)) , (6) 



X{x) = ^{x)U = -y=(i) 1 (x)+'il; 2 (x),'4) 1 (x)-'4)2{x)) (7) 

with 




U = -4s _ = ■ CO 



The discussion in this section proceeds via the following steps indicated by sub- 
sections. In subsection 2.1, the Grassmann integration is carried out exactly in two 
steps, first on the even and then on the odd sublattices. In this process, a procedure 
based on graphical elements (vertices) is established for describing the results. Com- 
binations of these graphical elements (vertex clusters) stand in correspondence to 
contributions to the partition function (sum) Z\. In subsection 2.2, rules for the cal- 
culation of the weight of the contribution of an arbitrary vertex cluster are discussed. 
In subsection 2.3, the fact that the weight of certain vertex clusters may vanish (due 
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to certain cancellations among different terms contributing to the partition function) 
is exploited to interpret the graphical rules established in subsection 2.1 in terms 
of a loop model picture. This finally allows to specify in subsection 2.4 the statisti- 
cal model equivalent to the one-flavour lattice Thirring model with Wilson fermions. 



2.1 The result of the Grassmann integration 

In the expression for the partition function (||) the integration on the even sublattice 
A e is now performed first. As each lattice site x supports 4 Grassmann variables 
terms which respect 



h + 2s x + 4t x = 4 (9) 

can give a non-vanishing contribution only (s x = 0, 1, 2 is the power of Mip(x)i(j(x) 
and t x = 0, 1 is the power of 2Gipi(x)ip2{x)ipi{x)ip2(x) in the expansion of the expo- 
nentials in eq. (§).)• Below we give the results for all local integrals related to a given 
lattice site x e A e which allow non-vanishing contributions. We have relied on a 



purpose- written Mathematica program [14| to do the necessary algebra. The results 



have been arranged in a way most suitable for the further discussion. The graphical 
rules below have to be interpreted as follows: the black dot denotes any point x on 
the even sublattice A e at which the Grassmann integration is performed, a dashed 
line means k\ = while thin and thick lines stand for ki = 1 and k\ = 2 respectively. 
A vertex is understood to be a lattice point with certain admissible values of the oc- 
cupation numbers k\ on the attached links. A given set of occupation numbers ki on 
the lattice A uniquely determines a certain diagram denoting a contribution to the 
partition function (As one can see from eq. (|i~0D, we consider the cases s x = 2, t x = 
and s x = 0, t x = 1 as unified into one.). Finally, the first coordinate component of 
x is understood to be the horizontal one while the second component runs vertically. 

There is one possible vertex with s x — 2 or t x — 1: 

Vertex 1 



I r 2 ( 1 

- + -- = J II d Mx) I - (M4j{x)^j{x 
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+ 2G if>i(x)ifa{x)il>i(x)ij>2(x)j 
M 2 -2G (10) 



All further vertices have s x = (vertices 2-20) or s x = 1 (vertices 21-26) and t x — 
(The results for the vertices 2-7 can also be found in the paper by Salmhofer ||, 
including the explicit steps of the derivation of the r.h.s. of eqs. flT2] ) and (|H|).). 



Vertex 2 



I 



Xi(x - e 2 ) xi{x + e 2 ) X2(x + e 2 ) \2{x - e 2 ) (11) 



Vertex 3 



^i(x - ei) i)i(x + ei) ?/> 2 (x + ei) ^ 2 (x - e x ) (12) 



Vertex 4 



= - ^2(2; + ei) % 2 (x - e 2 ) xi(x - e 2 ) Vi(a; + ei)/2 (13) 

Vertex 5 

— V— = - ipi(x - ei) Xi( x + e 2) Xz{ x + e z) ^{x- ei)/2 (14) 
1 
1 

Vertex 6 

^ 2 (x + ei) xi(^ + e 2 ) X2(a; + e 2 ) + ei)/2 (15) 

I 
I 

Vertex 7 



-r 



V>i(x - ei) % 2 (x - e 2 ) xi(a: - e 2 ) ^ 2 (x - ei)/2 (16) 
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Vertex 8 
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Vertex 9 



Vertex 10 



Vertex 11 



Vertex 12 



[V>i(x - ei) xi( x + e ?) Xi( x ~ e 2) i>2(x - ei) 

- - ei) x 2 (:r - e 2 ) x 2 (z + e 2 ) ^ 2 (^ - ei)]/2 (17) 



[V> 2 (z + ei) Xi( x + e 2 ) Xi(^ - e 2 ) ^i(a: + ei) 

- ^ 2 (x + ei) x 2 (x - e 2 ) x 2 (z + e 2 ) + ei)]/2 (18) 



[X 2 (x + e 2 ) + ei) ^i(^ - ei) + e 2 ) 

- X2{x + e 2 ) ^2(2: - e x ) 4) 2 {x + ei) Xi( x + e 2)]/2 (19) 



[Xi(x - e 2 ) Vi(a: + d) ipi(x - e x ) X2{x - e 2 ) 

- xi( x ~ e 2 ) ^2( x ~ ei) ^ 2 (x + ei) x 2 (:r - e 2 )]/2 (20) 



[V>i(z - ei) ^i(a; + ei) - e 2 ) ^2(2: - e x ) 

+ ^1(2; - ei) x 2 (x - e 2 ) ^2(2: + ei) ^2(2: - ei)]/v^ (21) 



Vertex 13 



— = [X2(x + e 2 ) ipi{x + ei) xi(x ~ e 2 ) + e 2 ) 

+ X 2 (a; + e 2 ) x 2 (a: - e 2 ) ^ 2 (x + d) + e 2 )]/>/2 (22) 

Vertex 14 

I 
I 

+ ^2(x + ei) + e 2 ) ^(x - ei) + ei)]/v^ (23) 

Vertex 15 

— = [Xi(^ - e 2 ) ^2(2: - ei) X2(a; + e 2 ) X2(a: - e 2 ) 

+ - e 2 ) xi(x + e 2 ) ip!(x - ei) x 2 (x - e 2 )]/v / 2 (24) 

Vertex 16 

[- ^i(x - ei) + ei) X2(a; + e 2 ) ^2(2: - e x ) 
+ ^i(x - ei) + e 2 ) ^ 2 (a; + e x ) ^ 2 (x - ei)]/v^ (25) 

- e 2 ) ^i(a; + ei) x 2 (a; + e 2 ) X2(x - e 2 ) 
- Xi(^ - e 2 ) + e 2 ) *p2(x + ei) x 2 (x - e 2 )]/v / 2 (26) 



Vertex 17 
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Vertex 18 



— = [~ X2{x + e 2 ) i) 2 {x - e x ) xi{x - e 2 ) Xi{x + e 2 ) 

+ X2(x + e 2 ) x 2 (x-e 2 ) i/j^x - e x ) Xi(x + e 2 )]/V2 (27) 

Vertex 19 

l 

— +— = [foix + ej i) 2 {x-e x ) Xi(x-e 2 ) $i(x + e{) 

- ^) 2 (x + e x ) X2{x - e 2 ) ipi(x - e x ) Tpx(x + ex)]/V2 (28) 

Vertex 20 

— " — = -[ Xi{x ~ e 2 ) i) X (x + ex) X2{x + e 2 ) i) 2 (x - ei) 

+ i)i(x - e x ) xi{x + e 2 ) xi{ x - e 2 ) $i(x + ei)/2 
+ ^)i(x - e x ) x 2 (x - e 2 ) X2^x + e 2 ) ^i(x + ei)/2 
+ ^) 2 (x + ei) xi{x + e 2 ) xi(x + e 2 ) $ 2 (x - ei)/2 
+ i) 2 {x + ei) x 2 {x - e 2 ) X2{x + e 2 ) $ 2 (x - e ± )/2 
+ ^i(x - ei) xi(x + e 2 ) ip 2 (x + ei) X2{x - e 2 )\ (29) 

Vertex 21 

--H--- = M [xi(x - e 2 ) xi(x + e 2 ) + X2(x + e 2 ) X2(x - e 2 )] (30) 
Vertex 22 

l 

—Jf— = M [ipi(x — ei) ipi(x + ei) + ip 2 (x + ei) $ 2 (x — ex)] (31) 
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Vertex 23 



= M [ X i(x - e 2 ) fa(x + ei) - ip 2 {x + ei) X2(x - e 2 )]/V2 (32) 



Vertex 24 



= M [^i(x- ei ) Xi(^ + e 2 ) - X2(x + e 2 ) ^ 2 (x-e 1 )]/V2 (33) 

l 
l 

Vertex 25 

= M [x 2 (^ + e 2 ) + + + ^i) + e 2 )]/>/2 (34) 



Vertex 26 

l 

= M [xi(x-e 2 ) Mx- ei ) + ^i(x-ei) x 2 (s - e 2 )]A/2 (35) 



In order to achieve further understanding it is useful to assign the remaining ■?/>, ip 
fields on the odd sublattice A G certain graphical symbols as follows (The full black 
dot denotes a point x on the even sublattice while the hollow circle stands for the 
point on the odd sublattice which the argument of the fields relates to.). 

O-— = 4> 2 (x — ex) (36) o— • = il>\{x — e-i) (37) 

= Mx + a) (38) — -O = Mx + e ± ) (39) 

I = X2(x-e 2 ) (40) J = xi(x-e 2 ) (41) 

| = X 2 (x + e 2 ) (42) J = XiOr + e 2 ) (43) 

As a rule, an arrow flowing out of a point of the odd sublattice (hollow circles) de- 
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notes a if?, x field variable while an arrow flowing into a point on the odd sublattice 
relates to a if?, x held variable. From this it is already clear that only those products 
of (four) Grassmann fields at any point of the odd sublattice allow non-vanishing 
contributions to the partition function Z\ which have two incoming and two outgo- 
ing arrows in their graphical symbols. We will refer to this fact as the vertex arrow 
rule. The symbols for the fields on the odd sublattice can now be used in a natural 
way to construct further graphical building blocks. As will become clear further 

D (for M = 1). 



below we choose signs in accordance with eqs. (|30|) 



o 



6 



o 



6 



+ 



o 



X2{x + e 2 ) X2(x-e 2 ) + Xi(x - e 2 ) Xi(x + e 2 ) 



(44) 



o 



-o 



o » « » o 



+ 



O ' • ' o 



ip 2 (x + ei) ip 2 (x - ex) + Vi(x - ei) ^i(x + ei) 



(45) 



1 

71 




i 

72 




Xi(x-e 2 ) if)i(x + ei)/V2 - ip 2 (x + d) x 2 {x - e 2 )/V2 (46) 




0i (x - ex) xi(x + e 2 )/V2 - X2^x + e 2 ) ip 2 (x - / V2 (47) 




M% + ei) Xi{x + e 2 )/V2 + X 2(x + e 2 ) ifj^x + e^/V^ (48) 
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1 1 

72 + 71 





= xi( x - e 2) ^{x - ei)/V2 + - e x ) X2(x - e 2 )/V2 (49) 

One immediately recognizes the significance of the numerical factors for each ver- 
tex. The rule simply is that each corner element (fIBD-(fi9|) is associated with a factor 
1 / v2. Using the above pictorial language for the remaining fields and their combi- 
nations on the odd sublattice we can gain now already some intuitive understanding 
for the structures appearing on the r.h.s. of eqs. ( |Tl]) -( p{f ) by interpreting them in 



terms of the graphical symbols introduced in eqs. (|36|) -(|43T), fl44]) -(|49D. We give some 
characteristic examples; similar relations can also be obtained for the omitted ver- 
tices. 



Vertex 3 



- o — • — o x o — • — o = o > • > o x o « • ' o (50) 



Vertex 6 




- 1 — o x 1 — o = - 

2 2 






(51) 



Since the vertices 8-20 correspond to sums of field products the following condensed 
notation is appropriate. Each of the pictures below however can also be expanded 
in terms of oriented line segments. 



Vertex 10 





(52) 



Vertex 14 




x 1— o (53) 
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Vertex 20 



+ 




o 



6 




X o 



+ 




X 



-o 



(54) 



It should be mentioned that most of the above relations do not generalize 

to the case when the Wilson parameter r is chosen not to be 1. However, interest- 
ingly enough eq. ([54]) still applies for r ^ 1 if for the graphical symbols (f4"4|)-([4"UD the 
expressions then obtained for the vertices 21-26 (eqs. @)-([35|)) are inserted (with 
M= 1). 



What remains to be done now is to discuss the result for the Grassmann inte- 
gration on the odd sublattice A D . For the moment we exclude vertex 20 (and its 
equivalent on the odd sublattice) from consideration as its complexity requires spe- 
cial attention. It will be discussed at a later stage, in subsection 2.2. Now, it is 
useful to introduce a graphical notation for the product of two field variables at a 
certain point x of the odd sublattice A Q . In view of the vertex arrow rule we only 
need products of two Grassmann fields where one of them carries a bar (i.e. oriented 
line segments in the pictorial language). We order the products always in such a 
way that for a given lattice point the out-arrow field variable stands first (For points 
on the odd sublattice this entails that the bar variable is the first factor in a product 
while for points on the even lattice it is always the second factor; see e.g. the order- 
ing applied in eqs. (|Tl"l)-(|3~5l).). We give just an example; all other combinations are 
defined in an analogous way (see eqs. (|36l)-(|4"3|) for the notation). 



ij)i(x) ipi{x) 



(55) 



It is also useful to establish a graphical symbol for the mass term combination of 
the fields. 



o 



M 



(56) 



Using these symbols we can now conveniently discuss the Grassmann integration on 
the odd sublattice A G . The graphical symbols below denote in an intuitive way the 
product of (four) Grassmann variables (and sums of them in eqs. (|63| ) - (|65|) ). The 
ordering is performed within pairs of fields which stand for an oriented line segment. 
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In each pair the bar variable stands first. As these products of two Grassmann 
variables are Grassmann even, no further ordering prescription is required. The 
black dot in the hollow circle of the point on the odd sublattice A Q (which is the 
central point in the diagrams below) indicates that the Grassmann integration has 
now been performed. In view of the vertex arrow rule one finds the following non- 
zero results (Dots which are close together denote one and the same lattice point). 
We give some characteristic examples only. 



tZ&ZZ = / II d M x ) d4> a (x) M x ) M x ) Mx) f 2 {x) = 1 (57) 

J a=l 




It is interesting to note that the results (|6l|), (|63|), (|65|) are symmetric under arrow 
exchange, while those for (|60|), (|62"D , (|64l) are antisymmetric. This consideration can 
also be extended to the vertices (]17|)-([2"8|), (|30|)-(|35|) on the even sublattice A e , when 
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the graphical symbols (]36|)-(f43|) are also used there. For the moment, we denote the 
different terms in (l30j)-(]35f) as subvertices. One finds that the subvertices 

on the even sublattice have exactly the same symmetry properties with respect to 
arrow reversal as their counterparts on the odd sublattice. 



2.2 The calculation of vertex cluster weights 

In order to gain further insight into the structures under study it is useful to con- 
sider in eqs. ([□])- (|2Hp, fl5U|)-(|55D the numerical factor in front of a given product 
of four Grassmann variables (where the sign is fixed by the ordering prescription 
used) as weight of a given subvertex (combination of oriented line segments) on the 
even sublattice A e . One then observes that the subvertices (products) in (pTT| ) - (|28|) , 
(|30|)-(p5D have exactly the same weights as their counterparts on the odd sublattice 
A Q (cf. eqs. (^7])-(^)). Consequently, from now on the distinction between even and 
odd sublattice vertices will no longer play any role (except in the discussion of vertex 
20 given further below). Furthermore, one can convince oneself that their weights 
can easily be given by means of a simple rule. Each oriented line segment can be 
associated with a certain weight (+1 for straight lines, ±l/>/2 for bent lines). The 
weight for any particular line segment can be read from the upper lines in eqs. (|44]) - 
(j49|). The weights in vertices 2-19 (eqs. (|TT|)-([2"8|)) can then be calculated by simply 
representing them by two oriented line segments (thick lines are understood as two 
thin lines with opposite orientations) and then multiplying the weights of these line 
segments. For example, this way one obtains from (fR]), fl47f ) the negative sign on 
the r.h.s. of (|13D , (|14D . To obtain the weights in (|30| ) - (|35|) one simply multiplies the 
line segment weights by another factor of M. The same rule applies to eqs. (|57|)-(l65|). 

We are now almost prepared to calculate the weight of any given vertex cluster 
(i.e. a connected diagram built of the vertices 2-26). One remaining subject to be 
discussed is whether our above ordering prescriptions for products of odd sublattice 
Grassmann fields related to the vertices on the even and odd sublattices are com- 
patible. As each cluster can be understood as consisting of a certain number of thin 
line loops we may study the ordering problem separately for a single loop (For this 
purpose we consider thick lines as two oppositely oriented thin lines.). With an eye 
to eq. (|5|) , one may now easily convince oneself that if for the even sublattice vertices 
the ordering prescription applied in eqs. (|TT|)- (|28|) , (|30|)-(|35D is used one need not 
to reorder the Grassmann variables related to the odd sublattice points in order 
to obey the ordering prescription used in eqs. (p^)-(|65|) except for one arbitrarily 
selected point in any single loop. Consequently, each thin line loop contributes an 
additional factor of -1, characteristic of any fermion loop in field theory, arising from 
ordering effects. Consequently, thick line loops do not contribute any minus sign 
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because they can be understood as two thin line loops. 

After the above discussion it is now straightforward to also analyze vertex clusters 
which contain the vertex 20 (and its analogue on the odd sublattice) we omitted so 
far. The remarkable eq. (]54]) immediately suggests how to proceed. Each vertex 20 
can be split into three different diagrams where the vertex is interpreted within the 
loop picture according to the graphical symbols on the r.h.s. of eq. ([54]), and then 
one may proceed exactly as for vertex clusters without vertex 20 (In total we have 
(3x4=) twelve, in part identical, products of Grassmann fields with different loop 
picture interpretations attached to them, i.e. four products for each of the three 
terms on the r.h.s. of eq. (0).). However, eq. (]54]) stands for a vertex 20 on the 
even sublattice A e and it remains to be seen whether an analogous interpretation is 
possible for the analogue of vertex 20 on the odd sublattice A Q . The results for it 
read as follows (We here apply the same conventions as in eqs. (p7])-(|5^|).). 




In the above equations for ordering purposes we have applied a certain interpretation 
in terms of oriented line segments. However, each of the above six pictures allows 
two different interpretations in terms of oriented line segments (leading in total 
again to twelve different graphical terms). Consequently, one may decide to apply 
one ordering in a certain fraction a G R of cases (i.e. for those cases shown in d66|). 
flBTp) and in another fraction (1 — a) the other, alternative prescription. The weight 
for each subvertex then has to be multiplied by a or (1 — a) respectively in order to 
obtain the correct final result for the partition function. The two vertices in eq. (|66| ) 
can both be interpreted in an alternative way (relative to the one shown in eq. (j66j) ) 
as two oriented bent line segments passing each other without intersection. The 
four vertices in eqs. (|6"T| ) each allow an interpretation as either two oriented bent 
line segments passing each other without intersection (as shown) or as two oriented 
straight line segments with an intersection. In order to make convenient contact with 
our above discussion for clusters without the vertex 20 (and its analogue on the odd 
sublattice), it now turns out that for eq. (^) the choice a = 1/2 is to be made while 
for eq. (|67j ) a = — 1 has to be selected^ If one does so the following simple rule for 

1 Clearly, the final result for the partition function does not depend on the choice of a; however 
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treating the vertex 20 (and its analogue on the odd sublattice) emerges. Replace in 
any vertex cluster under consideration each vertex 20 (and its analogue on the odd 
sublattice) by two line segments according to the graphical rule 




Each of the line segments can be oriented, consequently each of the three terms 
on the r.h.s. of eq. ( |68|) yields four different diagrams of two oriented line segments 
(twelve different terms in total). While clusters without the vertex 20 (and its ana- 
logue on the odd sublattice) contain self-avoiding (ki = 1) thin line loops only, the 
last term on the r.h.s. of eq. ( |5"%D now allows (ki = 1) thin line loops with intersec- 
tions also. The weight of any of the twelve possible diagrams containing now instead 
of the vertex 20 two oriented (hi = 1) thin line segments, passing each other without 
intersection or intersecting each other, can simply be calculated by taking separately 
the weights of each oriented line segment from eqs. (f44l) -(f49D and multiplying them. 
This concludes the discussion of the vertex 20. 

The rule for calculating the weight of a vertex cluster consequently is the follow- 
ing. Replace any vertex 20 (and its analogue on the odd sublattice) according to the 
rule (|68| ) by two (hi = 1) thin line segments (If P20 is the number of its occurrences 
in a given vertex cluster, this rule leads to S^ 20 different configurations. The latter 
term from now on always denotes the diagrams generated from a vertex cluster by 
application of the rule (|68D, or ([72]).). Provide any of the, say, m (ki = 1) thin line 
loops (without reinterpreting thick lines as two parallel thin lines) with an orienta- 
tion. This leads for a given configuration to 2 m different graphs (From now on the 
latter term denotes diagrams generated this way from configurations.). Represent 
the obtained graphs completely as a collection of, say, n oriented thin line loops 
(with thick lines interpreted as two oppositely oriented thin lines). In total, this 
leads for any original vertex cluster under consideration to different S^ 20 2 m graphs 
with a certain number of oriented thin line loops each (this number may differ from 
graph to graph) . Calculate the weight of each single oriented thin loop of length I by 
multiplying the weights of the / combinations of oriented line segments it is built of 
and attach the ordering factor -1. The weight of each graph is given by the product 
of the n oriented thin line loops it can be represented by. Any vertex cluster conse- 
quently has as weight the sum of the weights of the 3^ >2 ° 2 m graphs generated from it. 



the present choice admits a particularly simple way of describing contributions to the partition 
function. 
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In order to make further progress, first we have to study the weight of an arbitrary 
oriented thin line loop L. According to eqs. fl44|)-(49) the absolute value of the weight 



is given by the number of corners C(L) of the loop as well as the number of occur- 
rences of the vertices 21-26, denoted by Pm, and reads 2" 

-C(L)/2 

M Pm . According 

to the above considerations the sign of the loop weight is given 

where C_ (L) is the number of (oriented) right-up, up- right segments (cf . the second 

term on the r.h.s. of eqs. (fig), fl47f )). Now, it is useful to take into account the 

equation 

C_(L) = q(L) + 1 (mod 2) , (69) 
where q(L) is the number of self-intersections of the loop L. A proof of eq. flBTSp can 



be found in Appendix A. From eq. (|69|) immediately follows that the weight of a 
thin line loop under consideration is given byQ 

(-l)Q( L ) 2~ C ^/ 2 M Pm . (70) 

This formula immediately also applies if L is understood to be a graph of, say, n 
thin line loops (q(L) can still be understood as the number of intersections; however 
only self-intersections of loops are really important because the number of intersec- 
tions of different loops is always even.). This graph of loops can exhibit 2 m different 
combinations of orientations if the configuration (obtained from a vertex cluster un- 
der consideration by application of the rule (^)) it is generated from contains m 
(ki = 1) thin line loops. However, q(L) may differ for these cases. Consequently, 

cancellations may and will occur among the 3 P20 2 m graphs generated from a given 
vertex cluster. In the following subsection we continue with some observations in 
this respect. 



2.3 From oriented to coloured loops 

For the moment, in order to simplify the discussion let us only consider vertex clus- 
ters without any vertex 20. We show that any vertex cluster L under consideration 
containing a (ki = 1) thin line loop with an odd number of the vertices 8-19 has 
vanishing weight. Let Pt be the number of the vertices 8-19 in a certain (k[ = 1) thin 
line loop A in L. Now, generate from any configuration with Pt{A) = 1 (mod 2) 
according to the procedure explained above 2 m graphs of oriented thin line loops 
(with thick lines interpreted as two oppositely oriented thin lines). There are, say, 
Pti vertices 8-19 in the (ki = 1) thin line loop A whose attached thick lines start 

incidentally, Pm is always even. This can easily be seen for any given vertex cluster by simply 
dropping all thick lines from it, i.e. by replacing all vertices 8-19 by vertices 21-26. 
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and end at the loop A, consequently Pti = (mod 2). If these thick lines are 
interpreted as thin lines, the number of thin line intersections is zero (modulo 2). 
There are Pt2 vertices 8-19 in the (k\ = 1) thin line loop A whose attached thick 
lines start at the loop A and end at some other (ki = 1) thin line loop. If these thick 
lines are interpreted as thin lines, the number of thick lines starting from the loop 
A and containing one or zero (modulo 2) thin line intersections are Pt2o and Px2e 
respectively. It holds 

P T {A) = P T1 + P T2o + P T2e = 1 (mod 2) . (71) 

Consequently, either Pt2o or Px2e is odd (and the other is even, of course). If one now 
reverses the orientation of just the (ki = 1) thin line loop A (with all other orienta- 
tions fixed) the numbers Pt2o and Px2e interchange (and Aq = APt2o = 1 (mod 2) 
for the graphs with the opposite orientation of A). From eq. ( [70] ) one immediately 
recognizes that after summation over the two opposite orientations of the loop A 
the weight of any configuration containing A vanishes. 

We now draw some first conclusions from the insight just obtained. We demon- 
strate that any vertex cluster L can be built by the overlap of two sets of self-avoiding 
loops (Self- avoidance here means that also intersections of different loops within a 
given set are not allowed.). We distinguish the loops out of the two different sets 
by attaching a colour, say red and blue, to each sort of loop. We proceed as fol- 
lows. Divide any given (fcj = 1) thin line loop of a vertex cluster L into sections 
of alternating colour where the colour changes at any vertex 8-19 (This is possible, 
because any (ki = 1) thin line loop contributing to the partition function contains an 
even number of vertices 8-19, as we have shown in the preceding paragraph.). This 
procedure yields two different colourings for each (ki = 1) thin line loop in L (They 
just differ by an interchange between red and blue.). These two colourings stand 
in correspondence to the two different orientations of any (ki = 1) thin line loop in 
L. We can then interpret any thick line in L as two parallel thin lines of different 
colour. It is clear that the described procedure yields a covering of the given vertex 
cluster L by two sets of self- avoiding, differently coloured loops (at each vertex meet, 
at most, two links of the same colour only). Consequently, in view of eq. ([70]) , the 
weight of any vertex cluster (not containing any vertex 20) which is not zero must 
be positive. The graphs generated from vertex clusters with positive weight can be 
thought of as being constructed from the overlap of two differently coloured sets of 
self-avoiding loops (with monomer weight z = M and bending rigidity rj = l/y/2). 
Their weight is given by the product of the weights of the coloured loops they are 
built of. The number of different graphs corresponding to a vertex cluster is exactly 
given by the number of different coverings of it by means of two sets of differently 
coloured self-avoiding loops. 
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We would now like to extend the discussion to vertex clusters containing vertices 
of type no. 20. The best strategy to do so seems to be to transform any vertex cluster 
containing vertices of type no. 20 into "equivalent" vertex clusters without these 
(The problem we are interested in is primarily a topological one.). By "equivalent" 
we mean vertex clusters which have the same weight as the original one up to some 
positive constant which can easily be taken care of. The idea is to "stretch" any 
vertex 20 into two parts in such a way that in the "equivalent" cluster both parts 
are connected by a thick line. For any given column (row) of lattice points which 
contains a vertex 20 the "stretching" can be facilitated by inserting two additional 
lattice point columns (rows) to the left (top) and right (bottom) of the column (row) 
under consideration (The procedure is to be performed for each individual vertex 20 
separately). Pictorially, the procedure can be represented by the rule (for columns) 

— > J 1— + —I 1— • (72) 



Any other vertex 2-26 in the given column can also easily be "stretched" (shifted 
to the left or right) by appropriately inserting further horizontal thick or (ki = 1) 
thin lines (and filling the remaining sites emerging in the "stretching" procedure 
by monomers). Although one might be tempted to do so, the rule ([72]) should not 
be confused with the rule (^) which has a different algebraic meaning in terms 
of combinations of Grassmann variables. It is however a way of dealing with the 
vertex 20 which is alternative to QBlf). Let us be more specific to illustrate (|T2"D. 
For example, for a vertex 20 on the even sublattice A e one finds by means of the 
above procedure the expression (x denotes the coordinate of the vertex 20 on the 
"unstretched" lattice, which on the "stretched" lattice becomes the middle point of 
the introduced thick lines. On the "stretched" lattice the Grassmann integrations 
at the points x — e x , x, x + e x have been carried out.) 
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- i[)t(x - 2ei) x 2 (x -ex- e 2 ) X2(a; + ei + e 2 ) + 2e 2 )/2 



- 2ei) + ei + e 2 ) ^ 2 (a; + 2ei) X2(a; - ei - e 2 )/2 . (73) 



The first four terms correspond to the first term on the r.h.s. of (|72|) while the last 
four terms correspond to the second term on the r.h.s. of ([721). If one shifts back 
each argument in (|73|) by either plus or minus e\ one reobtains the expression on 
the r.h.s. of eq. (|29|) . Let us also discuss the rule (|72|) from the point of view of 
a vertex 20 on the odd sublattice A Q . Some of the corresponding results read as 
follows (Again, the Grassmann integration has been performed at the points x — ei, 
x, x + ei, where x denotes the coordinates of the vertex 20 on the odd sublattice.). 




In these equations we again apply the graphical conventions as explained previously 
(cf. eqs. (|D-([|9|), (H)). The equations ([F|), ([F|) relate to the first term on the 
r.h.s. of (|72|) , analogous results can be given for the second term. Eq. (|74l) corre- 
sponds to eq. (|66D (Remember that for the ordering choice selected in eq. (^) we 
have applied the weight a = 1/2, consequently we find —1/2 on the r.h.s. of eq. 
(|74D.). Eq. (p^) corresponds to the last two terms in eq. fl67|) (However, this time 
for the ordering choice selected in eq. fl67|) we have applied the weight a = 0. Con- 
sequently, eq. (^) is exclusively interpreted as intersection of oriented straight line 
segments. The corresponding reordering of the Grassmann variables introduces a 
relative minus sign. Therefore, we find 1/2 on the r.h.s. of eq. (|75|).). 

The rule ( |72|) now allows to efficiently handle the colouring problem for vertex 
clusters L containing vertex 20 and also to calculate their weights. For any given 
vertex cluster the rule ( |T2"D introduces 2 20 different configurations free of any ver- 
tex 20. As discussed above the summation over orientations of (ki = 1) thin line 
loops for these configurations (i.e. the summation over graphs) lets only those sur- 
vive which have an even number of vertices 8-19 in any (ki = 1) thin line loop. The 
sum of the weights of these configurations with non-vanishing (and consequently 
positive) weights yields the total weight of any vertex cluster containing vertices 
of type no. 20. For any graph generated from a configuration with nonvanishing 
weight the described colouring procedure can be carried out and then the rule (|72| ) 
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can be reversed without any problem (It depends on the shape of the configura- 
tion whether both terms on the r.h.s. of ([F^) lead to diagrams with non- vanishing 
weight.). Consequently, the above insight for the calculation of the weight of vertex 
clusters without any vertex 20 immediately carries over to any configuration gen- 
erated from an arbitrary vertex cluster. Finally, it is clear that any overlap of two 
sets of self-avoiding loops generates a vertex cluster constructed out of the vertices 
2-26. The two possible orientations of (ki = 1) thin line loops correspond to the 
colour exchange among the two sets (or, more precisely, to the exchange of the loop 
configurations among the two sets of different colour). 



2.4 The two-colour loop model 



We arrive at the following picture for the one-flavour Thirring model with Wilson 
fermions. Any diagram (and its weight) contributing to the partition function Z\ 
can be regarded as being generated by the overlap of two differently coloured sets 
of self-avoiding loops with monomer weight z = M and bending rigidity r\ — \j a/2. 
Only the monomer weight M 2 — 2G (eq. fllPl) ) of the Thirring model differs from 
the product of the monomer weights of the two differently coloured self-avoiding 
loop models M 2 . However, for free fermions (G = 0) there is complete agreement. 
Consequently, the structure of the statistical model equivalent to the Thirring model 
with Wilson fermions is very close to that for a free Dirac fermion. Therefore, for 
G = we find 



M, 



1 

71 



(76) 



where Z\[z, 77] is the partition function for the self-avoiding loop model with monomer 
weight z and bending rigidity r\ as defined in [T^j (see also further below). This re- 
sult is not astonishing because the self-avoiding loop model with bending rigidity 
i] = 1/a/2 is known |16[ (see also |T7|), fIBfl to be a free fermion model |19[, [^(J 
(cf. also sect. 4 of |15]). The model has second order phase transitions (within the 
Ising model universality class) at z = M = 2 (this corresponds to k — 1/4) and 
z = M = (where it is equivalent to a 6-vertex model). At z = M = the central 
charge c of the self-avoiding loop model with bending rigidity rj — 1/ v^2 is c = 1 
p2[ while, as has been argued in [15[, at z = M = 2 it should have c = 1/2 
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which is rather obvious from eq. ( f76|) because any free Dirac fermion has c = 1. 



Let us now specify in detail the statistical model equivalent to the one-flavour 
lattice Thirring model with Wilson fermions (and Wilson parameter r = 1). This 
statistical model is a two-color loop model with bending rigidity r] = 1/ \/2 and 



22 



UJ\ Cl>3 UJ4 CU5 Uq UJ-j Co>g 



Figure 1: Vertices of the (one-colour) self-avoiding loop model (The intersection 
vertex is omitted because its weight u 2 is zero.). Loop segments are denoted by 
continuous lines. 



monomer weight M 2 — 2G. The vertices and weights of this model can best be de- 
scribed as follows. Consider the (one-colour) self-avoiding loop model with monomer 
weight z = M and bending rigidity r] = l/y/2 (For the vertices see Fig. 1.). The 
weights are 

uj\ = z = M , (77) 
u 3 = U4 = 1 , (78) 

U 5 = = L0 7 = LUg = 7] = -j= . (79) 

v 2 

This model is a special 8-vertex model (with U2 = 0, i.e. a 7-vertex model) and, 
therefore, we have applied here the standard 8-vertex model notation (for a discus- 
sion of the model also see [^5p . The vertices of the two-colour loop model equivalent 
to the Thirring model can now be found by taking two differently coloured sets of 
the vertices shown in Fig. 1 and overlapping them in any possible way. This leads 
to 49 two-colour loop model vertices (because each one-colour loop model has seven 
vertices with non- vanishing weight). The weights of 48 out of the total 49 of these 
two-colour vertices are simply given by the product of the weights of the one-colour 
vertices they are built of. The only exception is the overlap of the two one-colour 
monomer vertices (which is the two-colour monomer vertex) which has the weight 
M 2 - 2G instead of M 2 (cf. eq. (gDD). 



3 Discussion and conclusions 

The value c = 1 of the central charge c of the one-flavour Thirring model, which is 



independent of the coupling G, has been discussed by various methods in f2"3] , f24 



p5fl . With the concept of universality in mind, this result is not surprising since a 
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variation of G only affects the monomer weight of the two-colour loop model equiv- 
alent to the Thirring model, consequently at least in a certain vicinity of G = the 
Thirring model should be expected to exhibit the same value of the central charge 
as free Dirac fermions. 

Let us also comment on the relation of the above two-colour loop model to 
other statistical systems. Somewhat similar two-colour loop models have recently 



been discussed in p6fl . Furthermore, the two-colour loop model equivalent to the 
Thirring model can also be viewed as a 4-state vertex model. To see this note that 
the one-colour self-avoiding loop model, whose vertices are shown in Fig. 1, is a 
2-state model. This leads, by construction, to the understanding of the Thirring 
model as a 4-state 49-vertex model. Q-state vertex models have been studied, e.g. 
in p7|-||36|. However, the particular 4-state models discussed in fl2~8| , f35fl , [|36|1 do 
not contain the case of the 4-state vertex model equivalent to the Thirring model 
which seems to not have been investigated in the published literature on the subject 
(except, in view of eq. (|76f ), the free fermion case, where one may use the results 
obtained in P^j). Finally, it is interesting to note, just as an aside, that also certain 



other models with four states have recently been discussed within the context of 



fermionic systems in [3S], |39 



The present investigation also immediately allows to draw conclusions about 
another lattice model, namely the one-flavour strong (infinite) coupling Schwinger 
model with an additional four-fermion interaction term (i.e. the Thirring- Schwinger 
model with some finite four-fermion coupling constant G as in the Thirring model 
itself). For G = Salmhofer has found the exact equivalence of this model to a 
(one-colour) self-avoiding loop model with monomer weight z = M 2 and bending 
rigidity i] = 1/2 ||. From our consideration of the Thirring model immediately 
follows, that for G ^ one only has to replace the monomer weight z = M 2 by 
z = M 2 — 2G in order to obtain the corresponding equivalence. Consequently, the 
model can be discussed along exactly the same lines as done for G = in |T5] . 



To conclude, in the present article we have found the exact equivalence of the 
one-flavour lattice Thirring model with Wilson fermions and Wilson parameter r = 1 
to a two-colour loop model which can also be viewed as a 4-state 49-vertex model on 
the square lattice. The choice r = 1 has proved crucial in simplifying the analysis. 
It should be emphasized that the two-colour loop model is a statistical model with 
positive weights which allows its numerical simulation by standard methods. The 
loop model picture even might prove advantageous because it allows the application 
of efficient cluster algorithms (e.g. see [4~0|, [^T|). The particular loop model equiva- 



lence found for free fermions (cf. eq. (|76|)) might be helpful in the future analysis of 
multi-flavour four-fermion models in 2D. Also, it might serve as a suggestive starting 
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point for the further analysis of Wilson fermions in higher dimensions. 
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Appendix A 



In this Appendix we demonstrate that the relation 

C_(L) = q(L) + 1 (mod 2) (A.l) 

holds, where q (L) is the number of self-intersections of a given thin line loop L and 
C_(L) is the number of (oriented) right-up, up-right segments (cf. the second term 
on the r.h.s. of eqs. (|46|), (|47|)). We only consider loops which do not occupy a link 
twice which, in principle, is possible due to links with ki = 2 (thick lines). One may 
always put such loops which visit at least one link twice on a finer lattice where 
each link can be occupied only once. The analysis demonstrating the validity of eq. 
( |A. 1|) consequently can also be extended to such cases. 

We start with q = 0. The proof will be based on an induction with respect 
to the area Q(L) enclosed by the loop L. If the area Q(L) of L is one (in lattice 
units), i.e. \L\ =4 (\L\ is the length of the loop), one may convince oneself that eq. 
( [A. 1| ) holds (C_ = 1 in this case). Consider now a loop L with Q(L) > 1. There 
exists at least one loop L' with Q(L') = Q(L) — 1 whose interior differs from L by 
just one lattice square (This selection of L' is important because otherwise certain 
transformations to be discussed below would possibly not apply). Consequently, we 
may perform a deformation of the loop L' into the loop L. We will demonstrate 
that in such a transformation AC_ = C_(L') — C_(L) = (mod 2) holds. In the 
transformation a lattice square is added to the interior of the loop L' which has at 
least one boundary link belonging to the loop L' . There are three different cases 
to be considered: the lattice square under consideration either has one ( (|A.2p ), two 



|)) or three ( ( [A . 4| ) ) boundary links belonging to the loop V (In the case of two 
links we do not need to consider the case where two opposite sites of the lattice 
square belong to the loop L' because this would be in contradiction with our choice 
for V '.). Below we depict some of the situations that may occur. The complete set 
is obtained by rotating the pictures in (|A.2|) -( [A!4D in all possible ways. 



n 



(A.2) 




j 



(A.3) 
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(A.4) 



One may now provide the loop segments in ( |A.2j )-( |A~4| ) with an orientation [Ll is 
denoted by a continuous line.). Then, for each possible case one may explicitly 
convince oneself that in the transition from the loop L' to the loop L AC_ G 
{ — 2,0,2}. According to our induction assumption eq. (|A.1|) holds for the loop L' 
(Vt(L') = Q(L) — 1), consequently it also holds for L. Therefore, eq. ( |A.1| ) is valid for 
any loop with q = 0. Now we are prepared to discuss the case q > (I am indebted 
to C. Krattenthaler for the argument fE?|].). We perform an induction with respect 
to q. The case q = has just been considered. Next, assume a given oriented loop 
L has q(L) > self-intersections. Select a certain intersection and apply to it a 
reinterpretation in terms of two oriented line segments passing each other without 
intersection. This leads to two different loops A and B with q(A) < q(L) and 
q(B) < q(L) intersections each (q(A) + q(B) = q(L) — 1). As one may convince 
oneself easily, in this operation in three of the four different cases the total number 
of (oriented) right-up, up-right segments does not change. C_ only goes up by 2 in 
the case of the following intersection diagram. 



(A.5) 



Consequently, within the considered operation C_ remains constant modulo 2. Now, 
according to our induction assumption eq. ( [A. 1| ) is valid for the loops A and B, 
therefore 

C_ = [q(A) + 1] + [q(B) + 1] ee q(L) + 1 (mod 2) . (A.6) 
This proves eq. ( |A.1|) . 
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